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We reconsider the excitation spectra of a dimerized and 
frustrated antiferromagnetic Heisenberg chain. This model 
is taken as the simpler example of compiting spontaneous 
and explicit dimerization relevant for Spin-Peierls compounds. 
The bosonized theory is a two frequency Sine-Gordon field 
theory. We analize the excitation spectrum by semiclassical 
methods. The elementary triplet excitation corresponds to an 
extended magnon whose radius diverge for vanishing dimer- 
ization. The internal oscilations of the magnon give rise to 
a series of excited state until another magnon is emited and 
a two magnon continuum is reached. We discuss, for weak 
dimerization, in which way the magnon forms as a result of a 
spinon-spinon interaction potential. 



One dimensional antiferromagnetic systems have at- 
tracted a great deal of interest over the last decades. 
They are typical examples where the low dimensional- 
ity enhances quantum fluctuations producing effects com- 
pletely different than the ones expected by classical theo- 
ries of magnetism. The interest in one-dimensional mag- 
netic system have recently been strongly renewed in view 
of the discovery of a new class , in general inorganic com- 
pounds, with definite one-dimensional character . De- 
tailed spectral characterizations are now avaible and the 
theory is pushed to make precise predictions on how dif- 
ferent microscopic interactions influence the spectral re- 
sponse of this system. We discuss in this work the case 
of a gapped antiferomagnetic chain and analyze how an 
explicit dimerization and a frustrated interaction com- 
pete between them and how the magnetic spectra evolve 
with the microscopic parameters. A similar subject has 
recently been analized by different numerical technicsEJ,B. 
An analytical approach for the Gase of weak dimeriza- 
tion limit has been also proposed! We will discuss the 
similarities and differences with our results. W&take the 
recently discovered inorganic material CuGeOsQ as a ref- 
erence even though our result could be more general. Our 
started model Hamiltonian reads: 

4 = U 1 + <5 (- 1 ) i ) s * • s *+i + aS * ■ ^+2} (!) 



where J is the nearest-neighbor (NN) exchange coupling, 
a the frustration parameter and S measures the amount 
of the explicit dimerization. There are two different mo- 
tivations to associate this model-Hamiltonian with the 
magnetic spectra of a Spin-Peierls compound. First, the 
temperature dependence of the magnetic susceptibility of 
CuGeOa could only be accounted for if an important next 
nearest-neighbor (NNN) is included in the ID Heisenberg 
model used to describe this material.tilj It is therefore 
natural to include an explicit dimerization on the ex- 
change interaction to describe the magnetic excitations 
in the low temperature Peierls phase and this has beea.in 
fact the approach used in various theoretical studieslD,El. 
Second, if the adiabatic hypothesis is left out, the system 
turn to be a coupled spin-phonon one. Even though the 
magnetic interactions have mainly one-dimensional char- 
acter they interact with the three dimensional phonons. 
The excitations of the low temperature phase are strongly 
affected by the interchain elastic coupling. In fact, topo- 
logical solitons (kink) which are genuine excitations of 
the isolated chain cannot longer exist as free excitations 
because they create a zone with opposite dimerization 
phase relative to the surrounding chains. A kink and an 
antikink are strongly interacting excitations. Their inter- 
actions were ap/jupunted in previous work by a linear con- 
fined potentialdu or as producing a kind of domain config- 
uration where the kink and antikink oscillate around an 
equilibrium distance!!!! As it has been recently remarked 
the model Hamiltonian (|l|) is the simplest example where 
this processes could be studied. In this context the first 
term of Eq. ([!]) account for the explicit dimerization im- 
posed by neighboring chains and the second term for the 
tendency to spontaneous dimerization as we will discuss 
in the following. 

For a — and strong dimerization (5) the ground state 
is a product of singlet over the strong bonds. The first 
excitations correspond to promote a singlet to a triplet 
and then delocalize it to build states of definite mo- 
mentum. They form a band of spin-1 magnon excita- 
tions. In the opposite limit (5 = 0) translational sym- 
metry is not explicitly broken. However for a greater 
than a critical value (a c ~ 0.23) the system sponta- 
neously dimerize and a gap in the spectra is open. At 
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a = 1/2 (Majumdar-Ghosh point) the exact ground state 
is the double degenerate product of singlets over the nn 
bonds. The excitations have bean variationally evalu- 
ated by Shastry and SutherlandB (SS). They are mas- 
sive S=l/2 spinons which correspond to an uncoupled 
spin separating two regions of singlet dimers. The two 
previously discussed cases represent the extremal situa- 
tion of correlation length equal to the lattice constant. 
As we discuss below the structure of the lower energy 
triplet excitations remains similar for other values of the 
correlation length in the limits previoulsy discussed. 

The main question we address in this work is how the 
spectra evolve from one limit to another i.e. how the 
massive spinons bound in a magnon. To address this 
question we analyze this system in the limit of small S and 
a slightly larger than the critical value. In this region the 
low energy spectrum could be-studied by bosonization. 
The bosonized Hamiltonian ialJ: 

f uB 2 u 
H bos = J dx{^-U 2 + — (Q^) 3 + gi sincj) 

+3 2 cos(20)} (2) 

g% ~ S and 32 ~ a — a c . The proportionality constants 
are of the order of the unity depending on the short range 
cutoff in the bosonization procedure. 8 = v2x for the 
isotropic model and u is the spin wave velocity. The 
spectrum of this double frequency sine-Gordon (DFSG) 
field theory is not known in general. We start by analyz- 
ing both limits of g\ = and 32 = where the theory 
reduces to a single Sine-Gordon and a lot of information 
is disponible. In particular the-aemiclassical method of 
Dashen, Hasslasher and Neveullj gives the exact mass 
spectrum of the particles for this model. 

For 32 = (we are indeed analyzing the whole zone 
32 < because this marginal interaction renormalize to 
zero in this case) the excitation spectra consists of a kink, 
an antikink and two breathers. In the original spin lan- 
guage, the kink excitation carries S z — 1, the antikink 
S z = — 1 and the breathers S z = 0. The lower energy 
breather is degenerate with the kink and antikink and 
these three excitations give rise to a triplet branch, the 
correlate of the previously discussed magnon band. For 
31 = the only one particle excitations of the resulting 
SG theory are a kink and an antikink, no breather is 
found. This kink carries spin one-half and is the analog 
of the massive spinons of the SS variational wave func- 
tion. However note that in this parameter regime where 
the continuum approximation is reliable the characteris- 
tic width of these excitations is large compared with the 
lattice space. 

We are now in position to discuss the excitations in the 
intermediate zone. We start by implementing a semiclas- 
sical calculation of the spectra of the theory (Q). Even 
though the method is not expected to give exact results 
for this non integrable model, it gives a valuable non per- 
turbative information which will allow us to interpolate 
between the desired limits. The time-independent equa- 



tion of motion corresponding to @ is: 

2Q2 ( 9 x^o) + 3i cos (j) - 2g 2 sin 20 o = 
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The lowest energy (homogeneous) configuration is <fio = 
— ^ (mod 2n) (it is not degenerate with O = § as it 
would if gi = 0). Solitons are x-dependent solutions of 
the equation of motion with finite energy with respect to 
this homogeneous solution. Owing to the fact that the 
system is Lorentz invariant, moving solitons are obtained 
from the static ones by a Lorentz transformation. To look 
for these static configurations we write a first integral of 
Eq. (§: 



(1 + B)e 



(d x <Po) 2 — B sin O — j cos 2</> = cte (4) 



where we have introducted the quantity B = 31/(432) as 
a characteristic parameter to analize the evolution of the 
spectrum from one limit to another. £ = u/A and A = 
fl\/ u (9i + 432) is the mass gap to create a particle (let it 
call a meson) above the homogeneous vacuum. £ is the 
correlation length of the theory and measure the width of 
the kinks. The boundary conditions for finite energy are 
lim x _ > ± 00 </> = -§ (mod 2n) and lim x ^ ±00 d x (f> = 0. 
They fix the cte = B + |. Equation (Q) could now be 
integrated. The solution center at the origin is: 

t(x) = sin(^o) = 1 — 2 cosh 2 [xo /£] sech[(x — aro)/£] 

sech[(x + ar )/£] (5) 

with 
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The quantity xo define a length scale which can be iden- 
tified as the 'radius' of this excitations. In the limit of 
B << 1 we have xq » 1 and t(x) could be approxi- 
mated by tanh(a; — xq) tanh(x + xq), i.e. a kink-antikink 
pair of the system without explicit dimerization. In the 
opposite limit xq = and (|J) becomes a kink solution 
of the static sine-Gordon equation corresponding to the 
dimerized chain. So that, t(x) inter-poles between this 
limiting cases. 

4>q (x) = arcsin[t(x)] winds a complete round clockwise 
or counter-clockwise. Therefore this excitation carries 
S z = ±1. The S z — component of this triplet exci- 
tation could be taken into account in our semiclassical 
approach by including the periodic time-dependent solu- 
tion of the field equation which cannot be obtained by a 
boosting of the static solution t(x). After quantization, 
this solution will give additional singlet states which are 
not taken into account in the quantum states generated 
by the quantization of t{x). For the single Sine Gordon 
equation these are the breather-like solutions. Further- 
more, DHN have shown that the first breather is nothing 
but a renormalized meson. In this context the n-excitcd 
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breather state is considered as bound state of n mesons. 
For our DFSG equation there are not analytic expres- 
sions for these breather solutions. We will consider in 
the following that a renormalized meson correspond to 
the S z = component of our triplet excitation. That is 
to say, for the SU(2) invariant model we should expect 
that after an appropriated resumation of the pcrturba- 
tive series the meson mass A will be equal to the soliton 
mass to be introduced below. 

To discuss the physical content of solution || we show 
in Fig. (|IJ) the unstagered part of the local magnetiza- 
tion < S* > {^d x cj> in the bosonic representation) for 
different values of the ratio of the parameters B corre- 
sponding to g\ smaller, equal or greater than gi- The 
correlation length(£) has been fixed to be ten times the 
lattice constant in this Figure. The figure have been ob- 
tained by freezing t(x). This corresponds to the physical 
situation where two nonmagnetic impurities cut the chain 
and fix the dimerization phase at the border. This sit- 
uation have been recently studied by different numerical 
technicsEl. Note the similarities of the behavior showed 
in Fig. 1 and the one of Fig. 3a of ref. |l| For the clean 
translational invariant system the local magnetization is 
x-independent as the result of the combination of this 
pattern center at the different site of the chain. 



semiclassical expansion this is a rigid configuration which 
translates as a whole. 

We now include gaussian fluctuation around this clas- 
sical solution. The eigenvalues problem of the fluctua- 
tion operator corresponds to a Schroedinger-like equation 
which reads as: 

\-dl + V(x)}ip = E<p 
V{x) = -^[Bt+{l-2t{xf)] (7) 

where we have made the substitution x — > ^ . The eigen- 

frequencies are given by to = A\/E. Equation ^ has 
a zero energy solution (with eigenfunction ip = d x (f>o) , 
this is connected with the translational mode arising in 
the broken of translational invariance of t(x). We re- 
call that this is the only bound state of the fluctuation 
operator for the case of the single sine-Gordon solitons. 
For the present case we find in addition another bound 
state which split from this zero mode for finite B. For 
increasing B the energy of this mode increases going to 
the continuum spectra of V(x) for B — > oo. In Fig. [2| we 
show the evolution of this eigenvalue with B. This figure 
has been obtained by numerical integration of Eq. (u\) . 
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FIG. 1. Uniform part of the local magnetization carries 
by the configuration defined by t(x) for three values of the 
parameter B 

Only for very low B there are two well separated peaks 
corresponding to two spin-1/2 excitations at a distance 
of the order Xq/£. For increasing B this two peaks are 
washed up and a single spin one excitation emerges at 
the chain center. All the intermediate spins in this con- 
figuration are excited from their value in the ground 
state. The excitation we are dealing up is therefore an 
extended spin-1 magnon which dissociate in two indepen- 
dent spinons for very weak dimerization and collapse in 
the magnon of the dimcrized chain in the opposite strong 
dimerization limit. Note that at this leading order of our 
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FIG. 2. (a) The evolution of the non-zero eigenfrequency of 
the fluctuation operator as obtained by numerical integration 
of (Q). (b) The number of magnon branches predicted by the 
semiclassical calculation. The insertion shows a typical low 
energy spectrum. 

The excitation spectra of the theory in the sector of 
Sz = ±1 is spanned by the following states: 

• The state of the quantum particle (of mass M) 
builded around t(x). We identify this as belonging 
to a one-magnon branch with dispersion E(p) = 
uy/p 2 + M 2 

• The excited states of the magnon of mass M* — 
M + nu>. There are as much as adittional triplet 
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branches as nu reaches the value A where the con- 
tinuum of the Schroedinger equation starts. This 
continuum corresponds to: 

• Labeling by q the continuum of level of Eq. ((?]), 

they are u>(q) = u^J t q 2 + J This is just the kinetic 

energy of a meson with momentum q. Therefore 
this state corresponds to the scattering of a meson 
in presence of the soliton. Moreover, when one of 
the continuum modes is excited once, we get a two 
particle meson-soliton state. As we have associated 
a meson with the elementary S z — excitation 
of the theory these excitations correspond to the 
S z = 1 component of the two magnon continuum 
of our original spin chain. 

The semiclassical calculation predicts the apparition of 
additional equi-spaciated magnon branches. The num- 
ber of these additional branches diverges as B goes to 
zero; this is nothing but the two spinons continuum of 
the undimerized chain. For decreasing dimerization or 
increasing frustration the number of the excited magnon 
states decrease until value of B. Beyond this critical value 
only two magnon branches are found. In Fig. || (b) 
we show the behavior previously discussed. The inser- 
tion gives our prediction for the low energy spectra for 
B greater than the critical value. Previous exact diago- 
nalization studies of the spectra of Hamiltonian [l] have 
shown the apparition of an additional triplet brancbLm 
confirming this picture. The increase of the number of 
the triplet branches with 5 has also recently seen in nu- 
merical studies of this systemtl. 

Now, we address the question of the contact of this 
semiclassical approach and previous work assuming a lin- 
ear confinement potential between the spinons. 

For B « 1 (weak dimerization limit) it is customary 
to think the problem as the oiiejjpf two interacting kinks 
as has been recently proposeduld. If they are well sep- 
arated their interaction is given by an attractive linear 
potential, the one arising on the term proportional to gi 
in Eq. (||). At low energy the spectra of the system cor- 
responds to the one of two particle interacting via this 
linear potential. A ladder of bound states is obtained by 
solution of an effective Schroedinger equation. At high 
enough energy these bound states are truncated because 
it becomes favorable to create a new kink-antikink pair 
than to excite one of them. This is essentially the sce- 
nario discussed in ref. ( ||) 

Let us define an effective kink-antikink potential by 
the following procedure: take Eq.(|J) as a suitable kink- 
antikink form with xq a variable distance, replace this 
form in Eq. (||) so the Xo-dependence of the total energy 
gives our kink-antikink interaction energy. The result is 
show in Fig. (||) for two small values of B where the kinks 
are expected to be well defined excitations. 

Two competing interactions could be clearly identify. 
An attractive linear interaction and a repulsive one acting 
at small distance. The repulsive interaction is originated 



on the finite width of the kinks. The value of xq given 
in eq. @ represent the minimuun of the total potential. 
For very small B the slope of the linear potential is small 
and this minimum is not well defined. The kinks and 
antikinks will be always far apart and they never seen 
the barrier. Therefore the linear approximation for the 
kink-antikink potential will give essentially the same re- 
sult as a potential like the one showed in the figure. The 
semiclassical approach gives only poor quantitative re- 
sults in this zone because it describes the problem as the 
one of an harmonic oscillators with very small frequency. 
However the general qualitative behavior of the spectra 
is well reproduced, i.e. the divergence of the number of 
excited states in the limit of no dimerization. 

For increasing B, a well defined minimum appears 
in the potential energy. Now, at low energy the kink- 
antikink oscillates around this minimum as bounded by 
a spring. Our system behaves as a vibrating molecule 
made by a kink and an antikink. The semiclassical ap- 
proximation previously discussed precisely describes the 
quantum state of this harmonic oscillator. 

By further increasing B the kinks loose their identity, 
t(x) becomes almost a kink solution of the SG theory with 
#2 = 0. The semiclassical result represents an internal 
oscilation mode in which the kink waveform undergoes 
a harmonically varying shape change localized about the 
kink center. 

Note that even for very small B the kinks are not well 
separated in the order of their width, thererofe we should 
associate the excitation in a broad range of the parameter 
as a unique magnon dispersed over a great number of 
spins. 
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FIG. 3. kink-antikink interaction potential as denned in 
the text 

Finally we analyze the experimental situation in 
CuGeOa. Neutron scattering measurements of the mag- 
netic spectratJ shows only one dispersive excitation be- 
fore the continuum. The excited states previously dis- 
cussed has not been found in these experiments. Two 
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possible reasons could be given. _ As it has recently 
show by numerical diagonalizationo the spectral weight 
of these states are very small. Moreover, in the original 
spin-phonoHjnodel the inclusion of the interchain elas- 
tic coupling^ implies that spins interact with transversal 
acoustic phonons. Therefore this excited state could de- 
cay in acoustic phonons so broadening the magnon peak. 
More experimental and theoretical work will be needed 
to elucidate this point. 
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